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Abstract: In this contribution we will discuss recent results concerning the intensity and the angular 
distribution of the gamma-ray and neutrino emissions as should be originated from the hadronic scatter- 
ing of cosmic rays (CR) with the interstellar medium (ISM). We assumed that CR sources are supernova 
remnants (SNR) and estimated the spatial distribution of primary nuclei by solving numerically the dif- 
fusion equation. For the ISM, we considered recent models for the 3D spatial distributions of molecular 
hydrogen. Respect to previous results, we find the secondary gamma-ray and neutrino emissions to be 
more peaked along the galactic equator and in the galactic centre which improves significantly the per- 
spectives of a positive detection. We compare our predictions with the experimental limits/observations 
by MILAGRO and TIBET (for the gamma-rays) and by AMANDA-II (for the neutrinos) and discuss the 
detection perspectives for a km3 neutrino telescope to be built in the North hemisphere. 
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Introduction 

Several orbital observatories (see [3] for a review), 
especially EGRET [16, 8], found that, at least up 
to 10 GeV, the Galaxy is pervaded by a 7-ray dif- 
fuse radiation. While a minor component of that 
emission is likely to be originated by unresolved 
point-like sources, the dominant contribution is ex- 
pected to come from the interaction of galactic cos- 
mic rays (CR) with the interstellar medium (ISM). 
Since the spectrum of galactic CR extends up to the 
EeV, the spectrum of 7-ray diffuse galactic emis- 
sion should continue well above the energy range 
probed by EGRET. That will be soon probed by 
GLAST [13] up to 300 GeV and by air shower ar- 
rays (ASA) (e.g. MILAGRO [19, 18] and TIBET 
[12]) above the TeV. 

Above the GeV, the main 7-ray emission processes 
are expected to be the decay of ir° produced by 
the scattering of CR nuclei onto the diffuse gas 
(hadronic emission) and the Inverse Compton (IC) 
emission of relativistic electron colliding onto the 
interstellar radiation field (leptonic emission). It is 
unknown, however, what are the relative contribu- 
tions of those two processes and how they change 
with the energy and the position in the sky (this 



is so called hadronic -leptonic degenerary). Sev- 
eral numerical simulations have been performed in 
order to interpret EGRET as well as forthcoming 
measurements at high energy (see e.g. [2, 1, 21] ). 
Generally, those simulations predict the hadronic 
emission to be dominating between 0. 1 GeV and 
few TeV, while between 1 and 100 TeV a compara- 
ble, or even larger IC contribution may be allowed. 

The 1-100 TeV energy range, on which we fo- 
cus here, is also interesting from the point of view 
of neutrino astrophysics. In that energy window 
neutrino telescopes (NTs) can look for up-going 
muon neutrino and reconstruct their arrival direc- 
tion with an angular resolution better than 1°. 
Since hadronic scattering give rise to 7-rays and 
neutrinos in a known ratio, the possible measure- 
ment of the neutrino emission from the Galactic 
Plane (GP) may allow to get rid of the hadronic- 
leptonic degeneracy. 

In this contribution we discuss the main results of 
a recent work were we modelled the 7-ray and 
neutrino diffuse emission of the Galaxy due to 
hadronic scattering [9]. Our work improves previ- 
ous analysis under several aspects which concern 
the distribution of CR sources; the way we treated 



CR diffusion by accounting for spatial variations 
of the diffusion coefficients; the distribution of the 
atomic and molecular hydrogen. 

The spatial structure of the ISM 

In order to assess the problem of the propagation of 
CRs and their interaction with the ISM we need the 
knowledge of three basic physical inputs, namely: 
the distribution of SuperNova Remnants (SNR) 
which we assume to trace that of CR sources; the 
properties of the Galactic Magnetic Field (GMF) 
in which the propagation occurs; the distribution 
of the diffuse gas providing the target for the pro- 
duction of 7-rays and neutrinos through hadronic 
interactions. In the following we assume cylindri- 
cal symmetry and adopt the Sun galactocentric dis- 
tance r Q = 8.5 kpc. 

The SNR distribution in Galaxy 

Several methods to determined the SNR distribu- 
tion in the Galaxy are discussed in the literature 
(see e.g. that based on the surface brightness - dis- 
tance relation [6]). Here we adopt a SNR distribu- 
tion a distribution as inferred from observations of 
related objects, such as pulsars or progenitor stars, 
as done e.g. in [11] which is less plagued from 
sistematics and agrees with that inferred from the 
distribution of radioactive nuclides like of 26 Al . A 
similar approach was followed in [22] where, how- 
ever, the contribution of type I-a SNR (which is 
dominating in the inner 1 kpc) was disregarded. 

Regular and random magnetic fields 

The Milky Way, as well as other spiral galaxies, 
is known to be permeated by large-scale, so called 
regular, magnetic fields as well as by a random, or 
turbulent, component. The orientation and strength 
of the regular field is measured mainly by means 
of Faraday Rotation Measurements (RMs) of po- 
larised radio sources. From those observations it 
is known that the regular field in the disk of the 
Galaxy is prevalently oriented along the GR Fol- 
lowing [14, 15] we adopt the following analytical 



distribution for the disk and the halo: 



S rcg (r, z) = B cxp<^ - 



where B = BJ*j*(r ,O) 



tb J 2cosh(z/z r ) 
(1) 

. , 2 //( I is the strength 
at the Sun circle. The parameters vb and z r are 
poorly known. However we found that our final 
results are practically independent on their choice. 
In the following we adopted tb = 8.5 kpc and 
z r = 1.5 kpc. 

More uncertain are the properties of the turbulent 
component of the GMF. Here we assume that it 
strength follows the behaviour 



B ran (r,z) = cr(r) B reg (r,0) 



1 



2 cosh(z/z t ) 



• (2) 



where a(r) parematrise the turbulence strength. 
Here we assume z t — 3 kpc. From polarimetric 
measurements and RMs is known GMF are chaotic 
on all scales below L max ~ 100 pc. The power 
spectrum of the those fluctuations is also poorly 
known. In [9] we considered both a Kolmogorov 
(B 2 {k) cc fc~ 5/3 ) and a Kraichnan (B 2 {k) oc 
fc -3 / 2 ) power spectra. 

The gas distribution 

The model which consider here is based on a suit- 
able combination of different analyses which have 
been separately performed for the disk and the 
galactic bulge. For the H2 and HI distributions 
in the bulge we use a detailed 3D model recently 
developed by Ferriere et al. [10] on the basis of 
several observations. For the molecular hydrogen 
in the disk we use the well known Bronfman's et 
al. model [4]. For the HI distribution in the disk, 
we adopt Wolfire et al. [23] 2-dimensional model. 
In the following we will assume that helium is dis- 
tributed in the same way of hydrogen nuclei. 

CR diffusion 

The ISM is a turbulent magneto-hydro-dynamic 
(MHD) environment. Since the Larmor radius of 
high energy nuclei is smaller than L max , the prop- 
agation of those particles takes place in the spatial 



diffusion regime. The diffusion equation describ- 
ing such a propagation is (see e.g. [20]) 

-V i (D ij (r,z)V j N(E,r,z)) = Q(E,r,z) (3) 

where N(E, r, z) is the differential CR density av- 
eraged over a scale larger than £ max , Q{E, r, z) is 
the CR source term and Dij(E, r, z) are the spa- 
tial components of the diffusion tensor. In the en- 
ergy range considered in our work energy loss/gain 
can be safely neglected. Since we assume cylin- 
drical symmetry the only physically relevant com- 
ponents of the diffusion coefficients are D± and 
Da, respectively the diffusion coefficient in the di- 
rection perpendicular to B rog and the antisymmet- 
ric (Hall) coefficient. We adopted expressions for 
those coefficients as derived by Montecarlo sim- 
ulation of charged particle propagation in turbu- 
lent magnetic fields [5]. Respect to other works, 
where only isotropic diffusion was considered and 
a a mean value of the diffusion coefficient was es- 
timated from the observed secondary /primary ratio 
of CR nuclear species (see e.g. [21]), our approach 
offers the advantage to provide the diffusion coef- 
ficients point-by-point at any energy. We solved 
the diffusion equation using the Crank-Nicholson 
method by imposing N(E) = at the edge of the 
MF turbulent halo and by requiring that it matches 
the observed CR spectra at the Earth position for 
most abundant nuclear species. 



Mapping the 7-ray and v emission 

Under the assumption that the primary CR spec- 
trum is a power-law and that the differential cross- 
section follows a scaling behaviour (which is well 
justified at the energies considered in this pa- 
per), the 7-ray (muon neutrino) emissivity due to 
hadronic scattering can be written as 



f N a pp Y 7 (a) (4) 
ds Ip(E p ; r, z) n H (r, z) 



dn 7 {v) {E; b,l) 
dE 



Here I p (E p ; r, z) is the proton CR differential 
flux at the position r, z as determined by solving 
the diffusion equation; a pp is the pp cross-section; 
Y 1 ~ 0.04 and Y Vfi+Dli ~ 0.01 are the 7-ray and 
muon neutrino yields respectively as obtained for a 



proton spectral slope a = 2.7 (see Fig. 7 in [9, 7] 
and ref.s therein); the factor Jm — 1.4 represents 
the contribution from all other nuclear species both 
in the CR and the ISM; s is the distance from the 
Earth; b and I are the galactic latitude and longi- 
tude. Here we assume that the slope of the CR in- 
jection spectrum is 2.2 and that the GMF turbulent 
spectrum has a Rraichnan spectrum. Fluxes ob- 
tained by using different spatial structure and spec- 
tra of the turbulent magnetic field differ by a factor 
2 at most. In 1 we compare our results with those 
obtained in [2]. 
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Figure 1: The 7-ray flux profiles along the GP 
(left panel) and along I — (right panel) for 
E > 1 TeV, averaged over 1° x 1° angular 
bins. The continuos (red), dot-dashed (blue) and 
dashed (green) curves correspond respectively to: 
our work, a model with the same gas density distri- 
bution but a uniform CR flux, taken to match direct 
observations at the Earth; the model considered in 
[2]. The corresponding neutrino flux can be ob- 
tained by dividing this diagram by 3. 1 . 



Discussion 

First of all, we compare our results with EGRET 
observations in the 4—10 GeV energy range [8]. 
That is possible since, already for < 10 TeV, nu- 
clei propagation takes place deep into the spatial 
diffusion regime so that the behaviour of the dif- 
fusion coefficients do not change going to lower 
energies and our results can be safely extrapo- 
lated. As we showed in [9], our predictions are in 
good agreement with EGRET measurements along 
the GR but a small deficit which can be easily 
explained in terms of IC. Then we can reliably 
compare our results with measurements performed 
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above the TeV with air shower arrays (see Tab. 1) 
and NTs. We found that with the exception of 
Cygnus (where one or more sources are likely to 
increase the local CR density) in all other regions 
we predict fluxes which are significantly below the 
experimental limits. 

Concerning neutrinos, the only available upper 
limit on the neutrino flux from the Galaxy has been 
obtained by the AMANDA-II experiment [17]. 
Being located at the South Pole, AMANDA can- 
not probe the emission from the GC. In the re- 
gion 33° < I < 213°, \b\ < 2°, and assuming 
a spectral index a = 2.7, their present constraint 
is (> 1 TcV) < 3.1 x 1CT 9 (cm 2 s sr)" 1 . 

According to our model the expected flux in the 



same region is 



*„„+*„(> 1 TcV) ~ 4.2 x 
10 11 (cm 2 s sr) 1 . That will be hardly detectable 
even by IceCube. Slightly more promising are the 
perspectives of a km 3 neutrino telescope to be built 
in the North hemisphere. In [9] we found, however, 
that the Km3Net project may be able to detect an 
excess in direction of the GC only if a significant 
CR over-density is present in that region. 
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4GeV 
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3 TeV 
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